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Abstract 

We report on the calculation of virtual processes contributing to the production of a Higgs boson 
and two jets in hadron-hadron collisions. The coupling of the Higgs boson to gluons, via a virtual 
loop of top quarks, is treated using an effective theory, valid in the large top quark mass limit. 
The calculation is performed by evaluating one-loop diagrams in the effective theory. The primary 
method of calculation is a numerical evaluation of the virtual amplitudes as a Laurent series in 
D — 4, where D is the dimensionality of space-time. For the cases H — > qqqq and H — ► qqq'q' we 
confirm the numerical results by an explicit analytic calculation. 
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FIG. 1: (a) Lowest order process for vector boson fusion, (b) Example of a diagram contributing 
to the gluon fusion process in association with two jets. 

I. INTRODUCTION 

In this paper we study the production of a standard model Higgs boson in association 
with two jets. This is one of the most promising discovery channels at the LHC especially 
for a Higgs boson with a mass in the range 110 GeV< Mh < 180 GeV. At Born level there 
are two classes of processes which contribute, as illustrated in Fig. [T] 

In Fig-EJa), the Higgs is produced via vector boson fusion, while in Fig.CJb) the coupling 
of the Higgs boson to gluons is mediated by a top quark loop. In the limit in which the 
mass of the top quark tends to infinity the coupling can be treated using an effective theory 
as described below. We shall refer to this process as the gluon fusion process. Notice that 
the external gluons lines in Fig. E^b) could as well be replaced by quarks. 

The final aim of this study is the calculation of the Higgs + 2-jet rate, at next-to-leading 
(NLO) order, where the Higgs is produced using the effective coupling to gluons, 



In Eq. G?" is the field strength of the gluon field and H is the Higgs-boson field. The 
effective coupling A is given by 



where g is the bare strong coupling and v is the vacuum expectation value parameter, 
v 2 = (GpVty -1 = (246 GeV) 2 . The finite 0(g 2 ) correction to the effective operator has 



\a{1 + A)HG%G a 



(1) 
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been calculated 

Ug 

167T 



The full NLO result will require the evaluation of the virtual corrections to the Higgs + 4 



ation of the tree graph rates 
and the calculation of a 



parton processes, which are the subject of this paper, the calcu 
from the Higgs + 5 partons amplitudes already given in refs. 
set of subtraction terms. The subtractions remove singularities present in the real emission 
diagrams in the regions of soft and collinear emission. After integration over the momentum 
of the un observed parton they are added back to the virtual emission diagrams and cancel 
the singularities in those virtual terms. 

We believe this calculation would be a useful addition to the literature for several reasons. 
First, the effective Lagrangian approach appears to be valid for light Higgs boson mass if 
the transverse momentum of the associated jets is less than the top quark mass |a E|- 
Second, this process constitutes a 'background' to the experimentally interesting vector 
boson fusion process, Fig^a). A complete NLO calculation will improve knowledge of 
this 'background' process. In addition, because the vector boson fusion process has a well 
determined normalization, it is one of the most accurate sources of information about the 
couplings of the Higgs boson at the LHC g]. An uncontrolled background from gluon fusion 
process could compromise that measurement. For a comprehensive review of standard model 
Higgs physics, see ref. 

Note that the process calculated in this paper is distinguished from the vector boson 
process, Fig E^a), by the presence of colored particles exchanged in the t-channel. The 
exchange of color charge generates extra jet activity in the central region, allowing discrim- 
ination against this process by a jet veto. Although the efficacy of such a veto will finally 
have to be determined by experiment, it will still be interesting to see how this works at the 
parton level with a full NLO calculation 1 . 

In the large top quark mass limit, virtual corrections have been considered in the effective 
theory by previous authors. Loop corrections to the process H — > gg are considered at one- 
loop level in ref. ]j and at two loop level in refs. wj. The results for the processes 



H — > ggg and H — > qqg are given in refs. jl2L Il3j . In the following we shall describe results 



1 To a limited extent this has been looked at in ref. 3]. However in a tree graph calculation one cannot 
look at the effects of finite jet size or of the central jet veto. 
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for the virtual corrections to the processes 



A) 


H - 


-> qqq'q' 


B) 


H - 


-> qqqq , 


C) 


H - 


, 


D) 


H - 


-> , 



(4) 
(5) 
(6) 
(7) 



using the effective theory, Eq. (JTJ. 



II. LOWEST ORDER PROCESS 

A. H — ► ggg'g' 

We first perform the calculation of the matrix element for the process involving two 
distinct flavors of massless quarks, q and q' , process A, 



q(k l )+ q(k 2 )+q'(k 3 ) + <j{fa) . 



(8) 



At Born level, only the diagram in Fig. |2Ja) contributes. The color expansion of the ampli- 
tude can be written as 



M^(fa, k 2 , fa, fa) 



S^s 13 5 ll 8 13 

<4 12 «2 U 



a(°)(l,2,3,4) 



(9) 



where ij denotes the color index of the jth quark and we have introduced the notation 

a (0) (l, 2, 3, 4) = a^(fa, fa; fa, h 2 ] fa, h 3 ; fa, fa) , (10) 

where fa and hi denote the momentum and the helicity of quark i. The result for the squared 
matrix element summed over the spins and colors of the final state quarks and antiquarks 
is then 



A (fa,fa,fa,fa) = ) ^M^ (fa, fa, fa, fa) 



[S 13 S 24 - S23S14) + S12S34 + (f!3 - s 2i) + (Sl4 - S 23 ) 



S 2 S 2 



2S 34 S12 



The number of colors, N c , enters as V = A 2 — 1, so, for the case of SU(3), we have that 
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(a) (b) 

FIG. 2: (a) Lowest order process for H — > qqq'q' ■ (b) Second diagram for identical quark process 
H -> qqqq. 

V = 8. The Lorentz invariants are defined as s^- = {hi + kj) 2 = 2ki ■ kj. The momentum of 
the Higgs can be eliminated in terms of the four massless momenta, pn = —k\ — k 2 ~ ^3 — k 4 , 
so that 



M H = S X2 + S13 + S14 + S 23 + S24 + S34 • 



(12) 



B. H 



qqqq 



In the case of massless quarks of identical flavor, process B, 

H -> g(fei) + q{k 2 ) + g(/c 3 ) + g(fc 4 ) 



(13) 



the Born amplitude squared is determined by the two diagrams, shown in Fig. 12(a) and 
(b), which differ by the exchange of the final state anti-quarks. The color expansion of the 
amplitude can be written as 



M*(k 1 ,k 2 ,k 3 ,h) 



£1£S £1£3 

'4 12 jy «2 U 



^(1,2,3,4)- 
M*(k u k 2 , k 3 , fc 4 ) - M^(k l} k A , k 3 , k 2 ) . 



12 % 



a (0) (l,4,3,2) 
(14) 



The result for the matrix element squared, summed over the spins and colors of the final 
state quarks and antiquarks is given by 

Boiki, k 2 , k 3 , k 4 ) = ^ \ M o (h, k 2 , k 3 , h)\ 2 

= Ao(h, k 2 , k 3 , k 4 ) + Ao(ki, fc 4 , h, k 2 ) + B' (ki, k 2 , k 3 , fc 4 ) , (15) 

where the interference term is defined as 

B' (h, k 2 , k 3 , h) = -2 ^ Re \Mf(k u k 2 , k 3 , k 4 )* M^{h, fc 4 , k 3 , k 2 )] . (16) 



The result for the lowest order interference term, B' , is given by, 



$13 ~ •524) 2 (-Sl2'S34 + S14S23 — S13S24) 



2(si 3 S 2 4 + S14S23 - S12S34XS12S34 + S13S24 - Si 4 S 2 3y 
1 



■Sl2'Sl4'S23' S 34 



(17) 



with C f = (Nl - 1)/(2N C ) = 4/3. 



C. H -> qqgg 

We now turn to process C, 

H -> q(h) + q(k 2 ) + g(k 3 ) + g(k A ) . (18) 
At lowest order the amplitude is given by, 

M C = (T a3 T a 4)n . 2C (0) (1)2) 3 ;4) + (T a 4T a3 )M2C (0) (1;2) 3 )4) ? (lg) 

where 03, 04 are the color indices of the gluons and i\, i 2 are the color indices of the quarks. 
As before we have introduced notation of the form 

cS 0) (l,2,3,4) =cf\k 1 ,h 1 ;k 2 ,h 2 ;k 3 ,e 3 ;h,e i ) , (20) 

where is the polarization vector of gluon i and c\ (1, 2, 3, 4) = cf\l, 2, 4, 3). Explicit forms 
for the three independent helicity amplitudes can be found, for example in refs. P, Q| • The 
former reference also contains explicit results for the amplitude squared. 

D. H -> gggg 

Lastly we consider the matrix element for the process D, 

H - g{h) + g{k 2 ) + g(k 3 ) + g(h) . (21) 
At lowest order the four gluon matrix element has the structure 

M D = tr(TMDTM2)T a -(3)T a -( 4 )) 4 0) (^(1),^(2),(x(3),(t(4)), (22) 



where the sum runs over the six non-cyclic permutations and we have introduced the notation 



2, 3, 4) = difa, e x \ fa, e 2 ; fa, e 3 ; fa, e 4 ) . (23) 
The partial amplitudes satisfy the relations |3, [r| 

4 0) (1,2,3,4) = df(4,l,2,3) cychcity , (24) 

4 0) (1,2, 3,4) = 4° } (4, 3, 2,1) reflection, (25) 

df\l, 2, 3, 4) + df ] (2, 1, 3, 4) + df } (2, 3, 1, 4) = dual Ward identity , (26) 

mplitu 



itudes. Explicit 
The former 



so that at Born level for fixed helicities there are only two independent am 
expressions for the helicity amplitudes can be found for example in refs 
reference also contains explicit results for the amplitude squared. Eqs. (|24j) and (|25p continue 
to be valid beyond leading order 

III. HIGHER ORDER PROCESSES 

In order to control the divergences which will occur at higher order we will continue the 
dimensionality of space-time, D = 4 — 2e. Within the context of dimensional regularization 
there remain choices of the dimensionality of internal and external gluons which are needed 
to completely specify the scheme. The most commonly adopted choices are the conventional 

;he 't Hooft-Veltman scheme, (HV) ^J, and the four- 



dimensional regularization, (CDR) 
dimensional helicity scheme, (FDH) (19|,|20|. In the CDR scheme one uniformly continues all 
momenta and polarization vectors to D dimensions. The HV scheme differs in the treatment 
of the external states, which remain four-dimensional. Finally in the FDH scheme all states 
are four-dimensional, and only the internal loop momenta are continued to D dimensions. 

Since we are interested in numerical evaluation, it is preferable to consider the external 
quarks and gluons in four dimensions, with two physical helicity states. We choose to work 



in the 't Hooft-Veltman scheme. The re 
schemes has been presented in refs. 121 



ationship of the CDR, HV and FDH regularization 



221 ]. It is therefore straightforward to translate our 
results to another scheme. The details of the translation between the HV and FDH schemes 
are provided in Section HVI 
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A. Distinct quarks 



At next-to-leading order in the perturbative expansion, 30 virtual diagrams contribute 
to the amplitude given in Eq. (jSJ). At one- loop level the amplitude can be decomposed into 
two independent color structures, 



M A (k u k 2 ,k 3 ,h) 



S^s i3 - —S il S i3 



i i; (l,2,3,4) + ^^a«(l,2,3,4). (27) 



The color sub-amplitude cup does not contribute at next-to-leading order because the inter- 
ference with the color structure of the Born amplitude vanishes. 

Before renormalization we find for the squared matrix element, summed over spin and 
colors of the final state 



A^hMMM) = J2(\ M o+ M i\ 2 ~\ M i\ 2 ) 



A (ki, k 2 , h, k 4 ) 



l + ^-Y A (k u k 2 ,k 3 ,h) 



+ A 



2 V g* 



X A (k u k 2 , k 3 , ki) + X A (k 3 , fc 4 , k u k 2 ) 



2 8tt 2 

+ X A (k 2 , h, k A , k 3 ) + X A (h, k 3 , k 2 , h)j + 0(e) . 
All ultraviolet and infrared singularities are in the functions Y{ki, k 2 , k 3 , k^) given by 



(28) 



Y A (k 1 ,k 2 ,k 3 ,h) = -N ( 



c r /i 



2c 



c 2 

e 2 



su) 6 + (s 23 ) 



+ 



1 cr/i 

,2e 



2c 



+ 



N c 
2 



3C> - b 



si 2 y e + (-s 34 )- £ - 2(-s 13 )- £ + 2(-s 14 )- £ + 2(-s 23 ) 

20 152 



2(-s 24 ) 



S12) + {-S 34 

T 2 \ , t Jme 



-n f + — N c - 16C f 



9 



9 



J- | S_ 1 (^134, S234, S34, Mfj) + Ls 

Ls?.7(s 12 3, s 134 ; S13, M 2 H ) + Ls^ e (si 24 , s 234 ; s 24 , M 2 H ) 



- (N c -—) hs^ e (s 12i ,s lu ;s 14 ,Ml) + Ls^ e (s 123 ,s 2 u;s 23 ,Ml) 



where 



and 



c r = (4tt) 



T(l + e)r 2 (l-e) 



(47T) 



r(i-2e) r(i- e ; 

11A C 2n/' 



+ 0(e 3 ), 



(29) 



(30) 



(31) 
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As usual rif is the number of light flavors and /x is the scale introduced to keep the coupling 
constant dimensionless in D dimensions. 

The finite function X A (fa, fa, fa, fa) is given by 

2 

X A (k 1 ,k 2 ,k 3 ,fa) = Ls_i(si2,si 3 ;si2 3 )-^- fi(fa,fa, fa, fa) 

2 

+ Ls_i(si 2 , s 23 ; s 123 )(N c - — ) fc 2 , & 3 , fa) 

+ + AUM^i) / 2 (A;i, fc 2 , fc 3 , fc 4 ) + iV c L (^^) / 3 (A;i, k 2 , fa, k 4 ) 
N c -s 12 -s 12 

+ — L (— ^) f 4 (ki, fa, k 3 , fa) + ( — + — )L (— ^) fc(fa, fa, fa, k 4 ) 
N c -si 2 2 7V C -si4 

- — L (— ^) / 5 (A;i, fc 2 , ^4, fa) + N c \n(-^-) f 6 (fa, fa, k 3 , k 4 ) 
N c -s 13 -s 12 

+ N c ln(^i^) f 7 (fa, fa, k 3 , fa) + N c H—) fs(h, k 2 , fa, fa) 

+ — ln(— ^) f 9 (fa, k 2 , k 3 , fa) + — ln(— ^) /i (A:i, k 2 , fa, fa) 
N c -s 12 N c -s 13 

^ $^2 3 1 

- — ln( — — ) /i (fc 2 , fa, fa, fa) + (N c + — ) fi 2 (fa, k 2 , k 3 , fa) . 
N c -s 23 N c 



+ 7rzr( ln (— — ) +M— -^)) (fu{fa,fa, fa, Jet) - fu {fa, fa, fa, fa] 
2JM r V — Si 3 — s 14 / V 



(32) 



The special functions coming from the loop integrals, L ,Li,Ls_i and Ls^™ 6 are given in 
Appendix El The explicit expression for the kinematic functions /$ are given in Appendix iBl 
We note that the line-reversal symmetry (1^2 and 3 «-» 4) and the renaming property 
(1^3 and 2 <-> 4) are manifest in Eq. (j28j) . 

The UV divergences are removed in the MS-scheme by adding a counterterm A ct given 

by 

A ct (fa, fa, fa, fa) = -2—b —^A (k u fa, fa, fa) . (33) 

Additionally, there is a finite contribution, A^, coming from the effective Lagrangian, 
Eq. (J^), which is 

A^ifa, fa, fa, fa) = 2A A (fa, fa, fa, fa) , (34) 
where A is given in Eq. ©. 
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B. Identical quarks 



In the case of identical quarks, 60 diagrams contribute the next-to-leading order process, 
Eq. (|13|) . Before renormalization we find for the squared amplitude, summed over colors 
and spins, 

B^fa, fa, fa, fa) = J2{\ M o+ M i\ 2 -\Mf\ 2 ) (35) 
= A^fa, k 2 , fa, fa) + Ai(kx, fa, fa, fa) + B[(fa, fa, fa, fa) , (36) 



with A\ given in (J28J). The result for the interference term can be written as, 
B[(fa,fa,fa,fa) = B' (ki, fa, fa, fa) 



l + ^Y B (fa,fa,fa,fa) 



+ A Z V 



2,r 9 



X B (fa, fa, fa, fa) + X B (fa, fa, fa, fa) + X B (fa, fa, fa, fa 



+ X B {fa, fa, fa, fa) + X B (fa, fa, fa, fa) + X B {fa, fa, fa, fa 



Bi 



+ X B (fa, fa, fa, fa) + X B (fa, fa, fa, k 
where the function Y B contains all divergent terms 



+ 0(e) 



(37) 



Y B (fa, fa, fa, fa) 



i(s 24 )-* + (s n )-*) 



+ 



+ 



c r /i 



2c 



AU 2 



-S12)" + (-S 34 )^ + (~S U )- e + (S 23 )- E - (-S24)" 



-S13J 



c r /i 



2c 



4e 



(-6C f + 2b ) 



[s 12 )- e + (s u ) e + (s 23 )~ e + (-S34)" 



20n f 80N C 8 
1 H H 

9 9 iV r 



1 

+ * 



LS 2 _T(S134, S234I S34, Mr) + Lsi7(Sl23, S234I S23, M H 



2me / 



+ Lsi7(si24, S134; su, M-„) + Lsi™ c (s 12 3, s 124 ; s 12 , MJ) 



LsiT(si23, S134; sis, M^)+ Lsi7(si 2 4, s 234 ; S 24 , 



(38) 
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The finite function X B is given by 

X B (fa, k 2 , fa, fa) = - Ls_i(si2, s 24 ; S124) gi{fa, fa, fa, fa) + -^j 



- Ls_i(si 2 , s 23 ; S123) 02 k 2 , fa, fa) 



N? 



+ Li (~~) ^(^1, ^2, fca, ^4) (l + j^) 

+ L (— ^) g 4 (fa, fa, fa, fa) + L (— ^) g 5 {fa, k 2 , fa, fa) — 



N 2 

c 

4- Id ( —J flr 6 (A;i, fc 2 , fo, + In ^ J^ ) ^(^1, fca, fa) 



S12 

"Sl23\ 



-Sl2 



+ gs(fa, fa, fa, fa) + 



(39) 



where the functions g i: are given in Appendix We note that the result in Eq. (J33J) is 
symmetric under the exchange of (1 3) or (2 <-> 4). 

The counterterm renormalizing the ultraviolet divergences in the case of identical quarks 
reads 

B ct (fa, fa, fa, fa) = -2—b —-zB (k 1 , fa, fa, fa) , (40) 

e l07T z 

while finite contribution coming from the effective Lagrangian is 



-B fin (/ci, fa, fa, fa) = 2A B (ki, fa, fa, fa) . 



(41) 



C. H -> g^c/ 



At one loop the full amplitude is calculated from 191 Feynman diagrams which can be 
decomposed into the three color-ordered sub-amplitudes, 

M C = (Tas T a 4)h .Jl) (lj 2) 3) 4) + ( T «4 T a3 )jii2C (l) (1? 2> g ? 4) + ^^^(D (1? 2> ^ 4) (42) 

Bose symmetry requires that (4(1, 2, 3, 4) = (1, 2, 4, 3). 

The divergent parts of these one-loop amplitudes are given by 



0^(1,2,3,4) 



4 1} (M,3,4) 



^ 167T 2 

— 1 H 

e e 



■^((- S 24)- 



cf(l,2,3,4) 



'-sis)- e + (-s 34 )- e + 



AU 2 



-S12J 



(43) 



16tt 2 



^4 0) (1, 2, 3, 4) ((- Sl4 )- e + (s 23 r ~ (sn)-' - (s 34 y 



+ ^4 0) (1 ; 2 ; 3 ; 4)( - . -(-s J2 r -(-s,i)- 



(44) 
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The interference between the Born and the NLO amplitude is given by 



2 Re(M 1 °M c 



VNc 
2 
V 



Re[(c« + c^cf + cDl + FReltfW + c^)1 , (45) 



2N r 



(i)^(o) 



,(0)> 



with Cj = Cj(l, 2, 3, 4). Counterterms, analogous to those in Eqs. (jHSllHIJ) need to be included 
to obtain the full renormalized result. 

Numerical results, which are given in the following section, were generated using an 
extension of the method su gges ted in ref. 12311 . Analytic expressions for the Feynman graphs 
are generated using Qgraf |24| and Form [25(. The scalar and tensor integrals appearing 
in the amplitudes are reduced numerically using the Davydychev reduction for the tensor 
integrals 0] and a recursive procedure similar to the one proposed in 0| to reduce all 
scalar integrals to a small number of analytically known basis integrals. These are then 
evaluated numerically as a Laurent series in the e parameter 2 . The key point of this method 
is that a record is kept of all previously computed integrals, so that each scalar integral is 
computed only once. The result of our procedure is a numerical expression for the scalar 
and tensor integrals component by component each of which has a Laurent expansion in e. 



This method will be described in detail in a later paper 2£ | . Numerica 



methods have also been described in refs 



n a later paper INumerica 

. 0,0,0,0, 22, HQS, 



or semi-numerical 
3- 



D- H -> gggg 



At NLO the amplitude for process Eq. (|21j) requires the calculation of 739 Feynman 
diagrams, which can expanded in nine color sub-amplitudes 

Mf = tr(TMi)TM2)TM3)TM4)) (a (1), a (2), a (3), a (A)) 

(7GS4/Z4 

+ — tr(T ai T aa ) tr(T a3 T a4 ) d£\l, 2, 3, 4) 

Ac 

+ — tr(T ai T a3 ) tr(T a2 T a4 ) d£\l,3,2,4) 



+ --tr(T ai T a4 ) tr(T a2 T a3 ) 4^(1, 4, 2, 3) . 



(46) 



2 The numerical Laurent expansion technique was first used in ref. |27| . In a more general analytic context 
it was used by many authors before. 
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If we discard diagrams with internal quark loops we have the decoupling identity 17 1 

4 1) (1,2,3,4)= d?(a(l),a(2),a(3),a(4)). (47) 

CT&S4/Z4 



However, at NLO the d 2 terms in Eq. (|4bj) do not receive contributions from internal fermion 
loops. This can be easily shown by explicitly examining the diagrams with internal fermionic 
bubbles, triangles, and boxes. The general expansion can thus be simplified as a consequence 
of Eq. (g7D so that 

Mf = ti(T a °wT a °WT a °wT a ^)d^\a(l),a(2),a(3),a(4)) 
— \tr(T ai T a ' 2 ) tr(T a3 T a4 ) + tr(T ai T a3 ) tr(T a2 T a4 ) 



N, 



+ tr(T ai T a4 ) tr(T a2 T as ) (1,2,3,4) . 

Using Eq. (|48|) it can be shown that the result for the matrix element squared is 

N 2 c (Nl-l) 



(48) 



|M ^ + MfT- \M[ 



16 



{\df\a(l),a(2),a(3),am 

0SS4/Z4 



+ 2 Re 



4 0) (a(l), a(2), (7(3), a(A)f d?(a(l), a(2), a(3), a(4)) 



(49) 



Counterterms, analogous to those in Eqs. (}3lfl need to be included to obtain the full 
renormalized result. 

Numerical results for this matrix element squared were generated using the method de- 
scribed above. The pole structure for the color sub-amplitude d\ has the simple form 



d? (1,2,3,4) 



16tt 2 



N c 

-4[(-S 12 r + (S23) 



-S34) + (-S14) 



4 0) (1,2,3,4). 

(50) 



IV. NUMERICAL RESULTS 

In this section we present numerical results for the Born amplitude squared and for its 
interference with the one-loop matrix element for the four processes of interest, A, B, C and 
D. We use the following arbitrarily chosen, momentum configuration, where a Higgs boson 
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of unit mass decays into four well separated partons, (E,p x ,p y ,p z ): 



Ph = 


(-1.00000000000, 


0.00000000000, 


0.00000000000, 


0.00000000000) , 


h = 


(+0.30674037867, 


-0.17738694693, 


-0.01664472021, 


-0.24969277974) , 


k 2 = 


(+0.34445032281, 


+0.14635282800, 


-0.10707762397, +0.29285022975) , 


k 3 = 


(+0.22091667641, 


+0.08911915938, 


+0.19733901856, 


+0.04380941793) , 


fc 4 = 


(+0.12789262211, 


-0.05808504045, 


-0.07361667438, 


-0.08696686795) . 



(51) 



For each process, {A, B, C, D}, we introduce the quantities 
1 



X B 



57T 



X (h, k 2 , k 3 , k 4 ) , 

X x {k u k 2 , h, fc 4 ) - X (k 1 , k 2 , k 3 , fc 4 ) 



with X = A, B, C, D , (52) 



g 6 A 2 L 

which are independent of the value of the coupling constant. Thus Xb is the matrix element 
squared evaluated using the Born amplitude. Xv,n an d Xv,a denote the contributions of 
the interference between the virtual amplitude and the lowest order, as calculated from the 
numerical and analytical formulas. The unrenormalized results are given in Table H] for the 
scale choice \i = M H and the momenta of Eq. (|51|). 

The explicit results show that far from exceptional momentum configurations, where 
divergent inverse Gram determinants are known to spoil the accuracy of the numerical 
procedure, a relative accuracy of O (10~ 13 ) can be achieved. For processes C and D, where 
a full analytical result is not available, we verified that the answer satisfies the Ward identities 
to a similar relative accuracy. For process D we checked numerically that for rif = 0, the 
color amplitudes satisfy the decoupling identity, Eq. (J4TJ) . Close to exceptional momentum 



28, 



configurations, it is still possible to use a numerical approach 

We have also checked numerically that our results satisfy the following relationship be- 
tween the HV and FDH regularization schemes, 

,2 



(1) FDH, 



a, L (1,2,3,4)-4 1)HV (1,2,3,4) 



9 



16tt 2 



N, 



3 -W)^ 2 ^ 



4 1)FDH (1,2,3,4)-4 1)HV (1,2,3,4) = 
c« FDH (l,2,3,4)-c« HV (l,2,3,4) 



16tt 2 



Nc 
6 



2N r 



cS 0) (l,2,3,4) 



4 1)FDH (1,2,3,4)-4 1)HV (1,2,3,4) = 
4 1)FDH (1,2,3,4)-4 1)HV (1,2,3,4) = 



(53) 
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cr/e 2 


cr/e 


1 


A B 








12.9162958212387 


Av,n 


-68.8869110466064 


-114.642248172523 


120.018444115429 


A V ,A 


-68.8869110466063 


-114.642248172519 


120.018444115458 


B B 








858.856417157052 


Bv,n 


-4580.56755817099 


-436.142317955660 


26470.9608978346 


B V ,A 


-4580.56755817094 


-436.142317955208 


26470.9608978350 


C B 








968.590160211857 


Cv,N 


-8394.44805516930 


-19808.0396331354 


-1287.90574949112 


Cv,A 


-8394.44805516942 


-19808.0396331363 


not known 


D B 








3576991.27960852 


Dv,n 


-4.29238953553022 TO 7 


-1.04436372655580 -10 s 


-6.79830911471604-10 7 


D V ,A 


-4.29238953553022- 10 7 


-1.04436372655580 -10 s 


not known 



TABLE I: Numerical results for the Born amplitude squared, (Xb), and the numerical and analytic 
one-loop corrections, (Xv,n and Xv,a), to the four processes A, B,C, D, Eqs. (J3HZI)- 



Applying the finite renormalization which compensates for the difference between the ultra- 
violet regularization in the two schemes j^], we recover the expected difference between the 
two schemes due to the differing infrared regularization, 

,2 



where 



a 



(1) FDH 

l 



(l,2,3,4)-ai 1)HV (l,2,3,4) = ^% a<°>(l, 2, 3, 4) 



cS 1)FDH (l,2,3,4)-cS 1)HV (l,2,3,4) 
4 1} FDH (1,2,3,4)-4 1} HV (1,2,3,4) 



4vr 2 

g 2 



!N7r ;a + 7 9 )cS 0) (l,2,3,4) 
^7,4 0) (1,2,3,4), 



lq = — and 7 9 e-. 



The other two relations in Eq. are unchanged. 



(54) 



(55) 
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V. OUTLOOK 



In this paper we presented results obtained using a general, semi-numerical calculation of 
one-loop corrections. In order to establish the feasibility of the semi-numerical method, we 
computed all the one-loop corrections to Higgs plus four parton processes using an effective 
Lagrangian. We presented explicit results for a specific, non-exceptional phase space point. 
For practical applications of this method, one has to be able to treat exceptional momentum 
configurations also. The method of this paper can be extended to the treat those regions. 
A detailed description of the algorithm is presented in a separate work 28^ . 

The results presented in this paper generate two separate lines of research. The first is 
clearly the completion of the calculation of the Higgs boson plus two jet process at next-to- 
leading order. As indicated in the text all of the needed elements are now in place. 

The second development is the extension of these methods to calculate other one-loop 
processes which currently lie beyond the range of analytic calculation. Examples of processes 
of current experimental interest are diboson plus one jet (Vi, V2,j), tri-boson production 
(^l) ^2, V3) and vector boson plus heavy quark pairs (VQQ). 
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APPENDIX A: INTEGRAL FUNCTIONS APPEARING IN AMPLITUDES 



The integral functions appearing in the virtual corrections are presented in this appendix. 



Following closely the notation of ref. 



3811 we define 



L„ W = ^M, L lW = i^fl±i. (Al) 
1 — r 1 — r 

The above functions have the property that they are finite as their denominators vanish. 
Furthermore we define 

Ls_ 1 (s ) t; m 2 ) = Li 2 (l - + Li 2 (l - -L) + In In - ^— , (A2) 

m z m z _ m z _ m z ^ 

16 



where the dilogarithm is defined as usual as 

Li 2 (x) = - f X dz ln{l - z \ (A3) 
Jo z 

The function Ls_i is simply related to the scalar box integral with one external mass eval- 
uated in six space-time dimensions, where it is infrared- and ultraviolet-finite. 

The 'easy' six-dimensional box function with two non-adjacent external masses, 1711,1713, 
is related to the function Ls 2 ™ 6 

Ls 2 _T (s, t; ml, ml) = -Li 2 (l - f ) - Li 2 (l - f ) - Li 2 (l - & ) - Li 2 (l - f) 

+L i 2 (l-^i)-Iln 2 (-f) . (A4) 

This function has the property that it vanishes as s + t — m\ — m\ — > 0. The analytic 
continuation of these integrals is obtained adding a small positive imaginary part to each 
invariant, Sij — > Sij + is. 
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APPENDIX B: FUNCTIONS FOR DISTINCT QUARKS 



The kinematic functions for the virtual corrections to H — > qqq'q' appearing in Eq. 
are given below: 



fi(k 1 ,k 2 ,k 3 ,k 4 ) = — ttt^ + 



•S12S34 3S13S24 — S23 + S14S23 — S14 — s|g Sl4 s 23 
2sj 3 S12S34 2si 2 sf 3 S34 



2 / \2 

S 24 lSi 3 S 2 4 — •S14S23J S24 . •S14S23 /-r>i\ 

- 7: I 1 1 2 2 \ UL ) 

2si 2 S 34 SJ 2 S 34 S13 S13 

, /, , , ,n Si 2 S34(Sl2S 3 4 + S23(>24 + 2S 23 - S U )) + s| 3 (s 2 4 + S U ) 2 

f 2 {h, k 2 , k 3 , h) = —3 (B2) 



f 3 (k 1: k 2 ,k 3 ,h) = + £23(524 + su) "^ ' ^ ' + "" Z4 ' ^" 23 (B3) 

ZS23 ^■ S 12 S 34 ^12 



2S12S34 



(B4) 



f (l 7 7 7 ^ S34 C \ S24 + 2S23 + 5Si4 

hih, k 2 , k 3 , k 4 ) = -2 s 2 3(s 2 4 + si4) —z 

S23 2sf 2 s 34 

4s 2 4 + 6s 2 3 — 3Si4 

2^ 

n n 7 7 7 \ s 13 2s 2 3 S 24 S23 S14S23 S13S23 S13S24 s 34 /q-v 

f 5 {h,k 2 ,k 3 ,h) = 1 2 1 1 ( B5 ) 

Sl4 ^24 S14 S34 S24S34 S24S34 S14S34 S14 
, /, , , , x S12S34 , 4S23S24 + 2Si 4 S 2 4 - 3Si 3 S 2 4 + 3S M S 2 3 

fe(k 1 ,k 2 ,k 3l h) = + 

2Si 3 S 2 3 2Si 2 S 34 

+ ^k^_fu + l (B6) 

-"I2S13.S34 S13 Z 

Mifex, fe, fc 8 , h) = ^ + Sl ff 12 ~ 

2sf 3 S34Si2 

, /, , , , x S14S23 - S13S24 m „x 

f 8 (ki,k 2 ,k 3 ,h) = (B8) 

^Si2S34 

h{Kk 2 MM) = ^i^_^_^|4^_ 2 £24 +2 £ 1 4_ 1 

S12S23S34 S13S23 S12S13S34 S 2 3 Sl3 

2Sn 4 + 2s 2 3S 2 4 + 5Si 4 S 2 4 - 5Si 3 S 2 4 + 5Si 4 S 2 3 



(B9) 
(BIO) 

M^htrt = ' 4{ 1' " MviS: ^ (bid 

S12S23S34 

fn(ki, k 2 , k 3 , k A ) = ^£Hi (B12) 

S 12-S34 

/l2 (^ A*, k 3 , h) = gl3(Sl3(gl4 " S24) + 2Sl4g23) + ^ ■ (B13) 

2Si 2 S 3 4 2S12 
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APPENDIX C: FUNCTIONS FOR IDENTICAL QUARKS 



The kinematic functions for the virtual corrections to H — > qqqq appearing in Eq. (|35|) 
are given below: 



S13-S24 0? 3 + S24) sf 3 + 2si 4 S 2 3 - 2si 3 S 2 4 + S 



2 

21 



01 A*, fes, **) = -1 - - + 

4S12S14S23S34 4si 2 s 34 

. gig ~ 2Si 3 S 2 4 + Sgj + 2g 12 g 3 4 
4S14S 23 

/, , , , \ 1 s 24 S 14 S 23 , 3Si 4 S23S 2 4 , S 13S 2 4 ($13 + ^34) 

g 2 [k u k 2 , k 3l fa) = -7--. o h o h 



4 4si 3 8si 2 Sf 3 S34 8S12S13S34 4S12S14S23S34 
2s^ 3 + 3 S14S23 — 3 S13S24 + 4^24 3 S12S24S34 S 12 S 34 



8Si2S 3 4 8S13S14S23 8s^ 3 Si4S23 



S14S23 + si 2 s 34 2sjg - 3 S13S24 + 4^24 + 3 S12S34 
8sf 2 Sx4 8sf 2 4si2 8S12S14 8S14S23 



'13 

S 23 S 24 S23 (S14 + S23 + S24) S34 (S 23 + S24) S34 ; , 



8s 13 8S14S23 

g 3 (fa, k 2 , fa, fa) = + r^r^_ ^ 1 _ j^l + vr^_^^ + _3^_ (C3) 



, , , x 6S 2 3-3S 2 4 . S 24 (-S23 + 4S 24 ) «34 . 5 S 24 S34 

g 4 (k l} k 2 ,k 3} k 4 ) = + — + - C4 

8S12 8S12S14 4S14 8S14S23 

S14 i S24 S24 , 3 S 23 ( — 2Si4 + S24) S23(Si4 + S24) 



g 5 h, k 2 , k 3 , k 4 4si2 + 8S12 4si 2 si 4 + 8si 2 si 4 4s£xS 34 

_l_ -^34 _|_ -S34 _ 3 S 2 4S34 _ S12-S34 (CB) 
4su 4s 23 8S14S23 4s M S23 

/, , , , x 5 s 2 4 3 S13S24 + 4^24 - 3 si 2 s 34 ,„ * 

9e (h,k 2 ,k 3 , fa ) = - - -. — + - (C6) 

8 4s u 8si 4 s 23 

/ 7 7 7 7 \ 7 S14 S 24 S? 4 S23 ^4 S12S34 
P7(«i,« 2 ,K3,«4) = -- + -: h- h 



'12°34 



4si 3 4su 4si2Si 3 S 34 4si 2 S 34 4si 3 S 23 4si 3 Si 4 S23 
2 

21 



Sl3 s 24 + 2^24 — S12S34 



8S14S23 

/ 7 7 ; 7 \ s 12 S14 S\ 2 S13 (S14 — 2^24) S13S24 

5f 8 (fci, fc 2 , fc 3 , fa) = — h — h — h 



32s M 32si2 32s 2 3 64si 2 s 2 3 32si 2 si 4 



(C7) 



S13S24 (S13 + S24) S23 (2Si4 + S24) S14 + S23 S13 (S12 + S23 — 2^24) 



64si2Si 4 s 2 3 64si 2 si4 32s 34 64si 4 s 34 

( s 12 — 2Si 3 + S14) S24 S13S24 (Sl3 + s 24) Si 3 S24 (^13 + s 24) 



64s 2 3S 3 4 64si2Si 4 S 34 64si 2 S 2 3S 3 4 

S13S24 (Sl3 + S24) , 0l3Sl4 + 2Si2 (S14 + S23) + S23S24) S 34 ,„„x 

+ — 1 — • {^°) 

04Si 4 S 2 3S 3 4 04Si 2 Si 4 S 2 3 



19 



[1] S. Dawson, Radiative Corrections To Higgs Boson Production, Nucl. Phys. B 359, 283 (1991). 
[2] A. Djouadi, M. Spira and P. M. Zerwas, Production of Higgs bosons in proton colliders: QCD 

corrections, Phys. Lett. B 264, 440 (1991). 
[3] V. Del Duca, A. Frizzo and F. Maltoni, Higgs boson production in association with three jets, 



JHEP 0405, 064 (2004) |arXiv:hep-p h/0404013 . 

[4] L. J. Dixon, E. W. N. Glover and V. V. Khoze, MHV rules for Higgs plus multi-gluon ampli- 
tudes, JHEP 0412, 015 (2004) |arXiv:hep-tn7041~ 1092 . 

[5] S. D. Badger, E. W. N. Glover and V. V. Khoze, MHV rules for Higgs plus multi-parton 
amplitudes, JHEP 0503, 023 (2005) |arXiv:hep-th/0412275| . 



[6] V. Del Duca, W. Kilgore, C. Oleari, C. Schmidt and D. Zeppenfeld, H + 2 jets via gluon 
fusion, Phys. Rev. Lett. 87, 122001 (2001) |arXiv:hep-ph/0105129| . 

[7] V. Del Duca, W. Kilgore, C. Oleari, C. Schmidt and D. Zeppenfeld, Gluon-fusion contributions 
to H + 2 jet production, Nucl. Phys. B 616, 367 (2001) |arXiv:hep-ph/010 8030 . 

[8] D. Zeppenfeld, R. Kinnunen, A. Nikitenko and E. Richter-Was, Measuring Higgs boson cou- 
plings at the LHC, Phys. Rev. D 62 (2000) 013009 |arXiv:hep-ph/0002036| . 

[9] A. Djouadi, The Anatomy of Electro-Weak Symmetry Breaking. I: The Higgs boson in the 
Standard Model, |arXiv:hep-ph/0503172j 



[10] R. V. Harlander, Virtual corrections to g g — > H to two loops in the heavy top limit, Phys. 

Lett. B 492, 74 (2000) |arXiv:hep-ph/0007289| . 
[11] C. Anastasiou and K. Melnikov, Higgs boson production at hadron colliders in NNLO QCD, 

Nucl. Phys. B 646, 220 (2002) |arXiv:hep-ph/0207004| . 
[12] C. R. Schmidt, H —> g g g (g q anti-q) at two loops in the large-M(t) limit, Phys. Lett. B 413, 

391 (1997) |arXiv:hep-ph/9707448| . 
[13] V. Ravindran, J. Smith and W. L. Van Neerven, Next-to-leading order QCD corrections to 

differential distributions of Higgs boson production in hadron hadron collisions, Nucl. Phys. B 

634, 247 (2002) |arXiv:hep-ph/0201114| . 



[14] R. P. Kauffman, S. V. Desai and D. Risal, Production of a Higgs boson plus two jets in 
hadronic collisions, Phys. Rev. D 55, 4005 (1997) [Erratum-ibid. D 58, 119901 (1998)] 
|arXiv:hep-ph/961054l1 . 



20 



[15] M. L. Mangano, S. J. Parke and Z. Xu, Duality And Multi - Gluon Scattering, Nucl. Phys. B 
298, 653 (1988). 

[16] F. A. Berends and W. T. Giele, Recursive Calculations For Processes With N Gluons, Nucl. 

Phys. B 306, 759 (1988). 
[17] Z. Bern and D. A. Kosower, Color Decomposition Of One Loop Amplitudes In Gauge Theories, 

Nucl. Phys. B 362, 389 (1991). 
[18] G. 't Hooft and M. J. G. Veltman, Regularization And Renormalization Of Gauge Fields, 

Nucl. Phys. B 44, 189 (1972). 
[19] Z. Bern, A. De Freitas, L. J. Dixon and H. L. Wong, Super symmetric regularization, two-loop 

QCD amplitudes and coupling shifts, Phys. Rev. D 66, 085002 (2002) |arXiv:hep-p h/02022"7l] . 
[20] Z. Bern and D. A. Kosower, The Computation of loop amplitudes in gauge theories, Nucl. 

Phys. B 379, 451 (1992). 
[21] Z. Kunszt, A. Signer and Z. Trocsanyi, One loop helicity amplitudes for all 2 — > 2 pro- 
cesses in QCD and N=l supersymmetric Yang-Mills theory, Nucl. Phys. B 411, 397 (1994) 

|arXiv:hep-ph/9305239] . 

[22] S. Catani, M. H. Seymour and Z. Trocsanyi, Regularization scheme independence and unitarity 
in QCD cross sections, Phys. Rev. D 55, 6819 (1997) |arXiv:hep-ph/9610553| . 

[23] W. T. Giele and E. W. N. Glover, A calculational formalism for one-loop integrals, JHEP 
0404, 029 (2004) |arXiv:hep-ph/0402152| . 



[24] P. Nogueira, Automatic Feynman graph generation, J. Comput. Phys. 105, 279 (1993). 



[25] J. A. M. Vermaseren, New features of FORM, arXiv:math-ph/0010025 . 

[26] A. I. Davydychev, A Simple formula for reducing Feynman diagrams to scalar integrals, Phys. 
Lett. B 263, 107 (1991). 

[27] A. van Hameren, J. Vollinga and S. Weinzierl, Automated computation of one-loop integrals 
in massless theories, arXiv:hep-ph/0502165 . 



[28] R. K. Ellis, W. T. Giele and G. Zanderighi, Semi-numerical evaluation of one-loop corrections, 

Phys. Rev. D 73, 014027 (2006) [arXiv:hep-ph/0508308| . 
[29] T. Binoth, J. P. Guillet, G. Heinrich, E. Pilon and C. Schubert, An algebraic / numerical 

formalism for one-loop multi-leg amplitudes, arXiv:hep-ph/0504267. 
[30] T. Binoth, G. Heinrich and N. Kauer, A numerical evaluation of the scalar hexagon integral 
the physical region, Nucl. Phys. B 654, 277 (2003) |arXiv:hep-p"h7 0210023 . 



in 



21 



[31] A. Ferroglia, M. Passera, G. Passarino and S. Uccirati, All-purpose numerical evaluation of 
one-loop multi-leg Feynman diagrams, Nucl. Phys. B 650, 162 (2003) [arXiv:h ep-ph/0209219]. 
[32] A. Andonov et ai, SANCscope - v.O.p, |arXi"v:hep-ph/0411186 



[33] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato and Y. Shimizu, 
Grace at one-loop: Automatic calculation of 1-loop diagrams in the electroweak theory with 
gauge parameter independence checks, arXiv:hep-ph/0308080 

[34] E. de Doncker, Y. Shimizu, J. Fujimoto, F. Yuasa, K. Kaugars, L. Cucos and J. Van Voorst, 
Loop integration results using numerical extrapolation for a non-scalar integral, Nucl. Instrum. 
Meth. A 534, 269 (2004) |arXiv:hep-ph/0405098| . 

[35] Z. Nagy and D. E. Soper, General subtraction method for numerical calculation of one-loop 
QCD matrix elements, JHEP 0309, 055 (2003) |arXiv:hep -ph/0308127 . 

[36] F. del Aguila and R. Pittau, Recursive numerical calculus of one-loop tensor integrals, JHEP 
0407, 017 (2004) |arXiv:hep-ph/0404120| . 

[37] W. Giele, E. W. N. Glover and G. Zanderighi, Numerical evaluation of one-loop dia- 
grams near exceptional momentum configurations, Nucl. Phys. Proc. Suppl. 135, 275 (2004) 



arXiv:hep-ph/0407016 . 

[38] Z. Bern, L. J. Dixon and D. A. Kosower, One-loop amplitudes for e+ e- to four partons, Nucl. 
Phys. B 513, 3 (1998) |arXiv:hep-ph/9708239| . 



22 



